No Title   |POLVER05_1
f(Ca) = v*Cao-v*Ca+ra*V  #
f(Cb) = v*Cbo-v*Cb+rb*V  #
f(Cc) = -v*Cc+rc*V  #
f(Cd) = -v*Cd+rd*V  #
r1a = -k1a*Ca*Cb^2  #
r2c = -k2c*Ca^2*Cc^3  #
ra = r1a+r2a  #
rb = r1b  #
rc = r1c+r2c  #
rd = r2d  #
r1b = 2*r1a  #
r2a = 2/3*r2c  #
r1c = -r1a  #
r2d = -1/3*r2c  #
Cao = 2  #
Cbo = 2  #
k1a = 10  #
k2c = 15  #
V = 2500  #
v = 100  #
Scd = Cc/(Cd+0.001) #
Ca(0)=2
Cb(0)=1
Cc(0)=1
Cd(0)=0.1